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1  Problem  Statement 

We  studied  the  quantum  adiabatic  algorithm  for  combinatorial  optimization,  first  proposed 
by  [8].  The  algorithm  works  as  follows.  One  constructs  a  time-dependent  Hamiltonian  H(t ), 
0  <  t  <  T,  which  can  be  efficiently  implemented  on  a  quantum  computer.  H{i)  is  designed 
so  that,  at  time  0,  the  ground  state  of  H{ 0)  is  known  and  can  be  prepared  efficiently,  while 
at  time  T,  the  ground  state  of  H (T)  encodes  the  solution  of  some  combinatorial  optimization 
problem.  Furthermore,  H{t)  changes  “slowly”  over  time. 

The  idea  is  that,  if  we  prepare  the  system  in  the  ground  state  at  time  0,  and  allow  it  to 
evolve  until  time  T,  then  the  system  will  follow  the  ground  state  of  H(t ),  and  end  up  finding 
the  solution  to  our  optimization  problem.  This  will  work  provided  the  change  in  H{t )  is 
sufficiently  slow;  in  particular,  by  the  adiabatic  theorem,  it  suffices  to  have  T  poly(l/gmin), 
where  gm ;n  is  the  minimum  gap  between  the  two  lowest  eigenvalues  of  H[t )  (for  all  times  t). 
Thus,  the  spectral  gap  of  H(t)  essentially  determines  the  running  time  of  the  algorithm. 
(Other  factors,  such  as  the  norm  of  H(t)  and  its  derivatives,  also  play  a  role,  but  they  are 
usually  less  significant.) 

For  example,  suppose  we  want  to  minimize  the  function  /  on  the  Boolean  hypercube 
{0, 1}”.  The  adiabatic  algorithm  might  use  a  Hamiltonian  of  the  form: 

H{t)  -  (1  -  t/T)H0  +  (f/T)Hi 

n 

i= 1 

Hi  =  ^  /(zi,..Mz„)  |^i)(^i|  ®  \zn){zn\. 

zi,...,zn€{0,l} 

(i) 

Here,  the  system  consists  of  n  qubits,  ax  denotes  the  ax  Pauli  operator  acting  on  the  rtli 
qubit,  and  |0)  and  |1)  denote  the  up  and  down  az  eigenstates  for  a  single  qubit. 
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The  adiabatic  algorithm  was  originally  proposed  to  solve  problems  such  as  3-SAT,  but 
no  rigorous  analysis  of  the  running  time  is  known  in  that  case.  Rigorous  analyses  are  only 
available  for  very  simple  “toy  problems”  (see,  e.g.,  [22,  7,  18]).  These  give  some  insight 
into  the  algorithm — in  particular,  they  show  how  the  adiabatic  algorithm  can  succeed  by 
tunnelling  through  potential  barriers,  and  how  it  can  fail  by  getting  drawn  into  local  minima. 
But  it  remains  difficult  to  understand  the  performance  of  the  adiabatic  algorithm  on  real 
optimization  problems,  because  the  situation  there  is  so  much  more  complicated.  At  best, 
there  are  heuristic  arguments  for  those  cases  (see,  e.g.,  [20]). 

Our  goal  was  to  develop  new  techniques  for  rigorous  analysis  of  the  adiabatic  algorithm, 
with  the  hope  of  understanding  the  performance  of  the  algorithm  on  some  moderately  compli¬ 
cated,  quasi-realistic  problems.  A  rigorous  analysis  of  this  type  would  be  useful  in  confirming 
(or  refuting)  the  conclusions  obtained  from  heuristic  arguments.  We  remark  that  this  sort 
of  analysis  has  been  carried  out,  to  some  extent,  for  classical  algorithms  such  as  simulated 
annealing  (see,  e.g.,  [21,  9,  1,  5,  6,  4]). 

2  Summary  of  Results 

The  following  is  a  summary  of  our  main  results.  For  more  details,  see  the  annual  progress 
reports,  and  the  reprints  of  papers  that  were  submitted  to  ARO. 

2.1  The  Adiabatic  Algorithm 

We  obtained  a  simple  proof  of  the  adiabatic  theorem,  using  a  martingale  argument  [11]. 
Another  more  powerful  proof  was  discovered  independently  by  Ambainis  and  Regev  [2] .  These 
results  give  some  intuition  about  why  the  spectral  gap  is  significant:  a  small  spectral  gap 
corresponds  to  a  rapid  change  in  the  ground  state  of  the  Hamiltonian,  which  in  turn  implies 
a  longer  running  time,  to  ensure  that  the  system  remains  in  the  ground  state.  In  many 
cases,  the  system  seems  to  undergo  a  quantum  phase  transition,  where  there  is  an  “avoided 
crossing”  between  the  two  lowest  energy  eigenvalues,  coinciding  with  a  sudden  change  in  the 
ground  state. 

Unfortunately,  we  made  little  progress  in  developing  new  tools  for  rigorous  analysis  of  the 
problem.  The  main  difficulty  is  that,  generally  speaking,  the  ground  state  of  the  Hamiltonian 
does  not  have  a  simple,  explicit  description.  The  only  known  exceptions  are  when:  (1)  the 
Hamiltonian  has  a  symmetry  that  effectively  reduces  the  dimension  of  the  problem,  e.g.,  it 
is  symmetric  under  permutations  of  the  qubits  [8,  22];  (2)  the  Hamiltonian  has  a  special 
structure  that  can  be  diagonalized  exactly,  e.g.,  it  can  be  transformed  into  a  system  of 
noninteracting  fermions  [18];  or  (3)  the  ground  state  has  relatively  uncomplicated  correlations 
and  entanglement,  e.g.,  certain  1-D  spin  chains  that  can  be  analyzed  using  the  martingale 
method  [17]. 

For  comparison,  note  that  in  classical  simulated  annealing,  the  stationary  distribution  of 
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the  Metropolis  random  walk  does  have  a  simple  description,  even  in  a  complicated  system — 
it  is  the  Gibbs  distribution.  This  makes  it  possible  to  calculate  the  second  eigenvalue,  for 
instance  by  variational  methods,  and  thus  bound  the  spectral  gap. 

Perhaps  heuristic  arguments  and  numerical  methods,  developed  in  condensed  matter 
physics,  are  a  more  viable  approach  to  understanding  the  adiabatic  algorithm.  For  instance, 
Schutzhold  and  Schaller  have  an  interesting  interpretation  of  the  algorithm  in  terms  of  first- 
and  second-order  quantum  phase  transitions  [19].  Also,  it  may  be  possible  to  simulate  some 
instances  of  the  adiabatic  algorithm  using  non-exact  methods  such  as  matrix-product  states 
[3], 

2.2  The  Local  Consistency  Problem 

We  also  obtained  several  results  on  the  complexity  of  the  “Local  Consistency”  problem.  (This 
is  not  related  to  the  adiabatic  algorithm,  but  it  is  nonetheless  interesting.)  The  problem  is 
as  follows:  given  a  collection  of  density  matrices  pi, . . .  ,pm,  where  pt  describes  a  subset  of 
qubits  Ci  in  an  n-qubit  system,  with  G) |  <  0(1),  decide  whether  there  exists  an  n-qubit 
state  a  that  agrees  with  all  of  the  pi  on  the  subsets  Ci . 

In  a  certain  sense,  this  problem  is  complementary  to  the  Local  Hamiltonian  problem  (find¬ 
ing  the  ground  state  energy  of  a  local  Hamiltonian),  which  is  QMA-complete  [10].  We  showed 
that  Local  Consistency  is  QMA-complete,  via  a  novel  reduction  using  convex  optimization 
with  a  membership  oracle  [12].  We  also  showed  that  7V-representability,  an  important  problem 
in  quantum  chemistry,  is  QMA-complete  [16].  In  addition,  we  proved  a  structural  property 
of  Local  Consistency  (that  if  a  solution  exists,  it  can  be  chosen  to  be  a  Gibbs  state)  [13];  and 
we  identified  special  cases  of  Local  Consistency,  for  1-D  and  stoquastic  systems,  whose  com¬ 
plexity  appears  to  be  strictly  easier  than  QMA  [15].  These  results  are  described  in  greater 
detail  in  my  thesis  [14]. 

2.3  Publications  and  Demographic  Data 

The  following  information  is  for  the  period  from  8/1/2006  to  8/31/2007  (since  last  year’s 
progress  report).  Papers  during  this  reporting  period:  [16]  (published  in  journal),  [14]  (PhD 
dissertation).  Demographic  data:  1  full-time  graduate  student  supported  by  this  agreement. 
Inventions:  none.  Technology  transfer:  none.  Other  activities:  Yi-Kai  Liu  gave  a  talk 
at  RANDOM  2006;  visited  the  Institute  for  Quantum  Computation  (IQC)  at  Waterloo  in 
January  2007;  gave  an  invited  talk  at  QIP  2007;  and  attended  the  SQuInT  2007  workshop. 
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We  study  the  computational  complexity  of  the  A-representability  problem  in  quantum  chemistry. 
We  show  that  this  problem  is  QMA-complete,  which  is  the  quantum  generalization  of  NP-complete. 
Our  proof  uses  a  simple  mapping  from  spin  systems  to  fermionic  systems,  as  well  as  a  convex 
optimization  technique  that  reduces  the  problem  of  finding  ground  states  to  IV-representability. 
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The  central  theoretical  problem  in  the  field  of  many- 
body  strongly  correlated  quantum  systems  is  to  find  ef¬ 
ficient  ways  of  simulating  Schrodinger’s  equations:  it  is 
very  easy  to  write  down  those  equations,  but  notoriously 
difficult  to  solve  them  or  even  to  find  approximate  solu¬ 
tions.  The  main  difficulty  is  the  fact  that  the  dimension 
of  the  Hilbert  space  describing  a  system  of  N  quantum 
particles  scales  exponentially  in  N.  This  makes  a  direct 
numerical  simulation  intractable:  every  time  an  extra 
particle  is  added  to  the  system,  the  computational  re¬ 
sources  would  have  to  be  doubled. 

The  situation  is  not  hopeless,  however,  as  in  principle 
it  could  be  that  all  physical  wavefunctions,  i.e.,  the  ones 
that  are  realized  in  nature,  have  very  special  properties 
and  can  be  parameterized  in  an  efficient  way.  The  idea 
would  then  be  to  propose  a  variational  class  of  wavefunc¬ 
tions  that  capture  the  physics  of  the  systems  of  interest, 
and  then  do  an  optimization  over  this  restricted  class. 
This  approach  has  proven  to  be  very  successful,  as  wit¬ 
nessed  by  mean  field  theory  and  renormalization  group 
methods.  So  far,  an  efficient  variational  class  to  describe 
complex  wavefunctions  such  as  those  arising  in  quantum 
chemistry  has  not  been  found. 

One  of  the  basic  problems  in  quantum  chemistry  is  to 
find  the  ground  state  of  a  Hamiltonian  describing  the 
many-body  system  of  an  atom  or  molecule.  The  essen¬ 
tial  element  that  makes  typical  Hamiltonians  very  un¬ 
generic  is  the  fact  that  at  most  2-body  interactions  occur. 
This  implies  that  the  number  of  free  parameters  in  such 
Hamiltonians  scales  at  most  quadratically  in  the  number 
of  particles  or  modes,  and  hence  the  ground  states  of  all 
such  systems  form  a  small-dimensional  manifold. 

In  the  case  of  a  Hamiltonian  with  only  2-body  interac¬ 
tions,  the  energy  corresponding  to  a  wavefunction  is  com¬ 
pletely  determined  by  its  2-body  correlation  functions, 
and  as  a  consequence  the  ground  state  will  be  the  one 
with  extremal  2-body  reduced  density  operators.  This 
fact  was  realized  a  long  time  ago,  and  led  Coulson  jj]  to 
propose  the  following  problem:  given  a  set  of  N  quantum 
particles,  can  we  characterize  the  allowed  sets  of  2-body 


correlations  or  density  operators  between  all  pairs  of  par¬ 
ticles? 

If  the  particles  under  consideration  are  fermions,  this 
has  been  called  the  N -representability  problem  [2}.  Here, 
we  consider  the  reduced  density  operators  acting  on  pairs 
of  fermions,  and  we  want  to  decide  whether  they  are  con¬ 
sistent  with  some  global  state  over  N  fermions. 

An  efficient  solution  to  the  iV-representability  problem 
would  be  a  huge  breakthrough,  as  it  would  (for  example) 
allow  us  to  calculate  the  binding  energies  of  all  molecules. 
Therefore,  a  very  large  effort  has  been  devoted  to  solving 
this  problem  and  there  exists  a  large  literature  on  this 
subject  (see  e.g.  0,0,0]). 

Here  we  will  show  that  the  iV-representability  problem 
is  intractable,  as  it  is  QMA-complete  and  hence  NP-hard. 
The  complexity  class  QMA  (Quantum  Merlin- Arthur)  is 
the  natural  generalization  of  the  class  NP  (nondetermin- 
istic  polynomial  time)  to  the  setting  of  quantum  comput¬ 
ing.  Colloquially,  a  problem  is  in  QMA  if  there  exists  an 
efficient  quantum  algorithm  that,  when  given  a  possible 
solution  to  the  problem,  can  verify  whether  it  is  correct; 
here  the  “solution”  may  be  a  quantum  state  on  polyno- 
mially  many  qubits  |23|  ■  A  problem  is  QMA- hard  if  it  is 
at  least  as  hard  as  any  other  problem  in  QMA;  that  is, 
given  an  efficient  algorithm  for  this  problem,  one  could 
solve  every  other  problem  in  QMA  efficiently  j2f[.  We 
say  that  a  problem  is  QMA-complete  if  it  is  in  QMA  and 
it  is  also  QMA-hard. 

In  a  seminal  work,  Kitaev  0  proved  that  the  Local 
Hamiltonian  problem — determining  the  ground  state  en¬ 
ergy  of  a  spin  Hamiltonian  that  is  a  sum  of  5-body  terms 
(on  n  qubits),  with  accuracy  ±e  where  e  is  inverse  poly¬ 
nomial  in  n — is  QMA-complete.  In  fact,  it  was  later 
shown  that  this  problem  remains  QMA-complete  when 
restricted  to  2-body  interactions  (7,  3,  and  even  in  the 
case  of  geometrically  local  interactions  ,9j. 

Because  the  Hamiltonians  under  consideration  are  lo¬ 
cal,  the  corresponding  ground  states  have  extremal  local 
properties.  The  dual  problem  of  determining  the  ground 
state  energy  of  a  local  Hamiltonian  is  to  decide  whether 
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a  given  set  of  local  density  operators  can  be  realized  as 
the  reduced  density  operators  of  the  same  global  state. 
Checking  the  consistency  of  such  a  set  of  local  density  op¬ 
erators  pW),  i.e. ,  checking  whether  there  exists  a  global 
state  o  compatible  with  those  local  density  operators, 
is  a  QMA-complete  problem  itself  E3-  In  the  present 
paper,  we  will  use  very  similar  techniques  to  prove  that 
TV-representability,  which  is  the  fermionic  version  of  that 
problem,  is  also  QMA-complete. 

To  sketch  the  main  ideas  of  the  proof,  we  will  first  con¬ 
sider  the  classical  marginal  problem:  suppose  we  have  N 
random  variables  that  can  take  the  values  ±1,  according 
to  some  joint  probability  distribution,  and  we  define  the 
2-variable  marginal  distributions 

p(ai,a,j)  =  ^2  p(a1,a2,...,aN). 

Given  a  set  of  marginal  probability  distributions,  we 
would  like  to  check  whether  they  are  consistent,  i.e., 
whether  there  exists  a  global  probability  distribution 
whose  marginals  are  equal  to  the  ones  we  were  given. 

Suppose  that  there  exists  an  efficient  algorithm  (whose 
running  time  is  polynomial  in  N)  to  solve  this  problem. 
Then  it  would  be  possible  to  identify  ground  states  of 
all  Ising  spin  glasses.  The  strategy  would  be  as  follows: 
since  the  set  of  consistent  marginals  is  convex  and  the 
energy  is  a  linear  function  of  the  marginals,  the  problem 
amounts  to  minimizing  a  linear  function  subject  to  a  set 
of  convex  constraints  that  can  be  checked  in  polynomial 
time.  This  can  be  done  in  polynomial  time  using  the 
ellipsoid  method  mt  It  therefore  follows  that  an  efficient 
algorithm  for  checking  consistency  allows  to  find  ground 
states  of  Ising  spin  glasses,  a  problem  that  is  known  to  be 
NP-hard  0.  Hence  the  problem  of  determining  whether 
a  set  of  binary  marginals  is  consistent  is  itself  NP-hard. 

Let  us  now  formulate  the  TV-representability  problem 
in  the  context  of  quantum  chemistry.  Electrons  and  nu¬ 
clei  tend  to  arrange  themselves  such  as  to  minimize  their 
energy,  and  the  binding  energy  of  a  molecule  can  be  de¬ 
termined  by  calculating  the  minimal  energy  of  the  corre¬ 
sponding  Hamiltonian.  In  practice,  the  nuclei  are  fairly 
well  localized  and  can  be  treated  as  classical  degrees  of 
freedom,  and  the  wavefunction  of  the  N  electrons  can  be 
approximated  as  a  linear  combination  of  tensor  products 
of  the  d  single-particle  modes  in  the  system  (those  form 
the  basis  set). 

As  electrons  are  fermionic,  the  complete  wavefunction 
must  be  antisymmetric,  and  this  is  most  easily  taken  into 
account  by  working  in  the  formalism  of  second  quantiza¬ 
tion: 

l 

*1,  'd  =  0 

U  +  ...+!«,  =  AT 

Here  aj  is  the  creation  operator  for  the  j’th  mode,  and 
|fi)  represents  the  vacuum  state  without  fermions.  The 


creation  and  annihilation  operators  obey  the  following 
anticommutation  relations: 

{uj ,  &j }  —  0  —  {cq ,  ( ij}  {cq ,  u—  )  —  2 Sij . 

Note  that  we  restrict  ourselves  to  the  subspace  of  states 
with  exactly  N  fermions,  d  denotes  the  number  of  modes, 
which  is  typically  much  larger  than  N.  The  number  of 
degrees  of  freedom  is  (^) ,  which  grows  exponentially  in 
N  when  d  >  cN  for  some  constant  c  >  1. 

In  the  case  of  quantum  chemistry,  the  Hamiltonian 
typically  contains  only  one-  and  two-body  interactions 
between  all  modes,  and  so  it  can  be  written  as  a  linear 
combination  of  terms  of  the  form  ajaj  and  a\a^aiak- 

The  2-fermion  reduced  density  matrix  (2-RDM)  is  cal¬ 
culated  by  tracing  out  all  but  two  of  the  fermions: 

p(2)  =  tr3j...iivp(Ar). 

where  p^N^  is  a  mixture  of  states  \tp)  with  exactly  N 
fermions.  In  the  language  of  second  quantization,  the 
matrix  elements  of  the  2-RDM  are  given  by: 

(2)  1  /  t  t  \ 

PijH  =  jv(jv_i  )K0I°J0*}- 

The  IV-representability  problem  (with  d  modes)  can 
now  be  stated  as  follows.  Consider  a  system  of  N 
fermions  and  d  modes,  d  <  poly(iV).  (For  purposes  of 
complexity,  we  consider  N  to  be  the  “size”  of  the  prob¬ 
lem.)  We  are  given  a  2- fermion  density  matrix  p,  of 
size  d<'d~1'>  x  ,  where  each  entry  is  specified  with 

poly(TV)  bits  of  precision.  In  addition,  we  are  given  a 
real  number  (3  >  1/ poly(iV),  specified  with  poly(TV)  bits 
of  precision.  The  problem  is  to  distinguish  between  the 
following  two  cases: 

•  There  exists  an  7V-fermion  state  a  such  that 
tr3,... ,jv(c)  =  p.  In  this  case,  answer  “YES.” 

•  For  all  N-fermion  states  a,  Ht^.  ..,n(ct)  -  p||i  >  /?■ 
In  this  case,  answer  “NO.” 

If  neither  of  these  cases  applies,  then  one  may  answer 
either  “YES”  or  “NO.”  Note  that  we  do  not  insist  on 
solving  the  problem  exactly;  we  allow  an  error  tolerance 
of  (3  >  l/poly(IV).  We  use  the  l\  matrix  norm  or  trace 
distance,  ||A||i  =  tr \A\,  to  measure  the  distance  between 
a  and  p. 

We  will  show  that  7V-representability  is  QMA- 
complete.  The  proof  consists  of  two  parts.  First,  we  show 
that  any  2-local  Hamiltonian  of  spins  can  be  simulated 
using  a  2-local  Hamiltonian  of  fermions  with  d  =  27V. 
Using  techniques  of  convex  programming,  we  show  that 
an  oracle  for  IV-representability  would  allow  us  to  esti¬ 
mate  the  ground  state  energies  of  2-local  Hamiltonians; 
thus,  IV-representability  is  QMA-hard. 

Second,  we  show  that  iV-representability  is  in  QMA; 
specifically,  we  construct  a  quantum  verifier  that  can 
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check  whether  a  2-particle  state  is  IV-representable,  given 
a  suitable  witness. 

Let  us  first  show  how  to  map  a  2-local  Hamiltonian, 
defined  on  a  system  of  N  qubits,  to  a  2-local  Hamiltonian 
on  fermions,  with  d  =  2 N  modes;  this  is  the  opposite  of 
what  has  been  done  in  0  The  idea  is  to  represent  each 
qubit  i  as  a  single  fermion  that  can  be  in  two  different 
modes  ai,bp,  so  each  IV-qubit  basis  state  corresponds  to 
the  following  IV-fermion  state: 

\zi)  ®  ®  M  ~  (4)1-*1  (b\yi  ■  ■  •  (ajv)1-^  (4)*"  |fi) . 

(1) 

Also,  all  the  relevant  single-qubit  operators  should  corre¬ 
spond  to  bilinear  functions  of  the  creation  and  annihila¬ 
tion  operators  (this  construction  guarantees  that  opera¬ 
tors  on  different  qubits  commute).  This  can  be  achieved 
by  the  following  mapping: 

erf  albi+blm,  erf  <->  %  (b\ai  -  a\b^j  ,  of  <->  1  -  2 b\bz. 

(2) 

It  can  easily  be  checked  that  these  fermionic  operators 
obey  the  Pauli  commutation  relations  within  the  sub¬ 
space  spanned  by  states  of  the  form  ©• 

Recall  that  the  qubit  Hamiltonian  can  be  written  as  a 
linear  combination  of  Pauli  operators.  So  the  fermionic 
Hamiltonian  can  now  be  created  by  rewriting  these  Pauli 
operators  with  respect  to  the  above  mapping  (note  that 
only  bilinear  terms  occur,  and  hence  operators  on  differ¬ 
ent  sites  commute) .  Notice  that  if  the  qubit  Hamiltonian 
only  contains  2-body  terms,  then  the  fermionic  Hamilto¬ 
nian  will  only  contain  terms  with  at  most  2  annihilation 
and  2  creation  operators. 

The  only  thing  that  is  left  to  do  is  to  guarantee  that 
there  is  exactly  one  fermion  on  every  site  i.  This  can 
be  achieved  by  adding  the  following  projectors  as  extra 
terms  in  the  fermionic  Hamiltonian: 

Pi  =  (2a j a*  -  l){2b\bi  -  1). 

All  the  Pi  are  biquadratic  and  commute  with  all  the  op¬ 
erators  introduced  in  ©,  and  hence  the  complete  Hamil¬ 
tonian  will  be  block  diagonal.  By  making  the  weights  of 
these  projectors  large  enough  (a  constant  times  the  norm 
of  the  Hamiltonian,  which  is  at  most  polynomial  in  IV), 
we  can  always  guarantee  that  the  ground  state  of  the  full 
Hamiltonian  will  have  exactly  one  fermion  per  site  1 2 fill . 

Let  us  now  assume  that  we  have  an  efficient  algorithm 
for  IV-representability.  We  claim  that  this  allows  us  to 
efficiently  determine  the  ground  state  energy  of  a  2-local 
Hamiltonian  on  qubits,  a  problem  which  is  known  to  be 
QMA-hard  0-  We  start  by  transforming  the  Hamilto¬ 
nian  L/qubit  on  N  qubits  with  2-particle  interactions  into 
a  Hamiltonian  f/fermi  on  d  =  2 N  fermionic  modes  with 
2-particle  interactions  (as  described  above).  The  prob¬ 
lem  of  finding  the  ground  state  energy  of  L/fermi  can  be 


expressed  as  a  convex  program,  which  can  be  solved  in 
polynomial  time  using  an  oracle  for  IV-representability 
(this  is  what  we  will  show  below).  Our  approach  is  similar 
to  0  ,  though  in  this  case  we  are  dealing  with  fermions 
rather  than  qubits. 

The  basic  idea  is  to  construct  a  convex  program 
that  finds  a  2-particle  density  matrix  p  which  is  N- 
representable,  and  which  minimizes  the  expectation  value 
of  .fffermi-  This  program  has  polynomially  many  vari¬ 
ables,  and  it  is  easy  to  see  that  the  set  of  IV-representable 
states  is  convex,  and  (iLfermi)  is  a  linear  function  of  p. 

There  are  two  minor  complications.  First,  Ufermi  de¬ 
scribes  a  system  with  2 TV  modes  and  an  arbitrary  number 
of  particles,  whereas  IV-representability  pertains  to  a  sys¬ 
tem  with  exactly  N  particles.  But  this  is  not  a  problem, 
because  the  interesting  behavior  in  Hfeimi  occurs  in  a  sub¬ 
space  where  all  the  states  have  exactly  N  particles.  So, 
in  our  convex  program,  we  only  need  to  consider  states 
p  that  have  exactly  N  particles. 

This  lets  us  simplify  7Vfermi  as  follows.  Since  we  are 
only  interested  in  how  it  acts  on  IV-particle  states,  we 
can  write  a\aj  =  jyzy<4(XX  a'ka,k)a,j ■  So  we  can  assume 
that  all  the  terms  in  Hfermi  are  of  the  form  afajajafc. 

Second,  convex  optimization  algorithms  usually  re¬ 
quire  that  the  set  K  of  feasible  solutions  be  full¬ 
dimensional,  i.e. ,  K  cannot  lie  in  a  lower-dimensional 
subspace.  So  we  have  to  represent  the  state  p  in  such 
a  way  that  there  are  no  redundant  variables. 

Let  S  be  a  complete  set  of  2-particle  observables  0, 
and  let  l  =  |«S|;  note  that  £  <  poly(d)  <  poly(IV).  We 
represent  p  in  terms  of  its  expectation  values  as  =  tr (Sp) 
for  all  S  £  S;  let  a  £  denote  the  vector  of  these 
expectation  values.  We  define  K  to  be  the  set  of  all  a 
such  that  the  corresponding  state  p  is  IV-representable. 
Note  that  the  iV-representability  oracle  allows  us  to  test 
whether  a  given  point  a  is  in  K. 

We  write  our  Hamiltonian  in  the  form  H  =  TS'-’ 

(plus  a  constant  term);  the  coefficients  ys  can  be  com¬ 
puted  using  the  Gram-Schmidt  procedure.  It  is  easy 
to  see  that  tr (Hp)  =  XXesTSa.S-  Our  convex  pro¬ 
gram  is  as  follows:  find  some  a  £  K  that  minimizes 

/(«)  =  Eses7s<AS. 

We  can  solve  this  convex  program  in  polynomial  time, 
using  the  shallow-cut  ellipsoid  algorithm  (see  0  ,  and 
references  cited  therein).  We  mention  a  few  technical 
details.  The  algorithm  requires  some  additional  informa¬ 
tion  about  K :  namely  a  guarantee  that  I\  is  contained 
in  a  ball  of  radius  R  centered  at  0,  and  K  contains  a  ball 
of  radius  r  centered  at  some  point  p.  (Also,  the  running 
time  of  the  algorithm  grows  polynomially  in  log {R/r).) 
In  our  case,  we  can  set  R  =  \fl,  and  r  =  1/ poly(£)  [27| . 

There  is  also  the  issue  of  numerical  precision.  Our  or¬ 
acle  for  IV-representability  has  limited  accuracy — it  may 
give  incorrect  answers  for  points  near  the  boundary  of 
K.  On  the  other  hand,  we  only  need  an  approximate 
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solution  to  the  convex  program,  in  order  to  solve  the  Lo¬ 
cal  Hamiltonian  problem.  It  turns  out  that  l/poly(AT) 
precision  is  sufficient.  The  ellipsoid  algorithm  still  works 
in  this  setting;  see  [3  for  details. 

Note  that,  in  place  of  the  ellipsoid  algorithm,  we  could 
have  used  a  different  algorithm  based  on  random  walks 

in  convex  bodies  13  ;  this  was  the  approach  used  in  E3- 

However,  it  is  not  clear  if  we  could  use  one  of  the  interior- 
point  methods  for  convex  optimization;  these  methods 
are  substantially  faster,  but  they  usually  require  an  ex¬ 
plicit  description  of  the  constraints,  not  just  a  member¬ 
ship  oracle. 

This  completes  the  proof  that  iV-representability  is 
QMA-hard.  As  a  corollary,  we  have  also  proven  that 
estimating  the  ground  state  energy  for  local  fermionic 
Hamiltonians  is  QMA-hard.  Note  that  the  use  of  the  el¬ 
lipsoid  algorithm  to  reduce  a  convex  optimization  prob¬ 
lem  to  a  convex  membership  problem  is  not  new;  see  13 
for  references  to  related  work. 

Next,  we  show  that  Wrepresentability  is  in  QMA. 
That  is,  we  construct  a  poly-time  quantum  verifier  V 
that  takes  two  inputs:  a  description  of  the  problem  (that 
is,  p  and  /?);  and  a  “witness”  r,  which  is  a  quantum  state 
on  polynomially  many  qubits.  The  verifier  V  should  have 
the  following  property:  if  p  is  ./V-representable,  there  ex¬ 
ists  a  witness  r  that  causes  V  to  output  “true”  with 
probability  >  pi',  if  p  is  not  Wrepresentable  (within  er¬ 
ror  tolerance  (3),  then  for  all  possible  states  r,  V  outputs 
“true”  with  probability  <  po;  and  pi  —  po  >  1/ poly(AT). 

The  idea  is  that,  when  p  is  ./V-representable,  the  correct 
witness  r  consists  of  (multiple  copies  of)  an  iV-fermion 
state  <7  that  satisfies  tr3i...]jv(cr)  =  p.  Then  the  verifier 
can  use  quantum  state  tomography  to  compare  a  and  p. 

We  represent  the  ./V-fermion  state  a  using  d  qubits,  via 
the  following  mapping: 

(a4)n...(c4)!d|Sl)  <->  |ii)  ®  ...  ®  | U). 

Call  the  resulting  qubit  state  a.  We  use  the  Jordan- 
Wigner  transform  to  map  the  fermionic  annihilation  op¬ 
erators  to  qubit  operators: 

Cli  <->  Ai  =  —(®k<i  Ofc)  ®  1 0)  ( 1 1 ^ . 

Thus,  an  observable  O  =  ajaja;afc  +  a\a\ajai  is  trans¬ 
formed  into  O  =  A\A^AiAk  +  A\A\AjAi:  which  is  a 
tensor  product  of  many  single-qubit  observables  and  one 
four-qubit  observable. 

We  claim  that  the  expectation  value  (O)  can  be  es¬ 
timated  efficiently.  Without  loss  of  generality,  assume 
that  the  eigenvalues  of  O  lie  in  the  interval  [0, 1].  Then 
there  is  an  efficiently-implementable  measurement  which 
outputs  “1”  with  probability  (O)  @] .  By  repeating  this 
measurement  on  multiple  copies  of  the  same  state,  we 
can  estimate  (O);  for  our  purposes,  it  is  enough  to  have 
polynomially  many  copies  of  the  state. 


We  now  describe  the  verifier  V .  The  witness  r  consists 
of  several  (i.e. ,  polynomially  many)  blocks,  where  each 
block  has  d  qubits,  supposedly  representing  one  copy  of 
the  state  d.  On  each  block,  V  measures  the  observable 
Y/k  1 1) (1  |fc,  and  if  the  outcome  does  not  equal  N,  V  out¬ 
puts  “false.”  This  projects  each  block  onto  the  space  of 
A-particle  states. 

Next,  V  performs  measurements  on  each  block,  to  es¬ 
timate  the  expectation  values  of  <r,  for  a  suitable  set  of 
observables.  Then  V  checks  whether  they  match  the  ex¬ 
pectation  values  of  p.  One  problem  arises:  the  prover 
could  try  to  cheat  by  entangling  the  different  blocks  of 
qubits.  One  can  show  that  this  does  not  fool  the  veri¬ 
fier,  using  a  Markov  argument,  as  was  done  in  13-  This 
suffices  to  show  that  ./V-representability  is  in  QMA,  and 
hence  finishes  the  proof. 

What  can  be  said  about  the  complexity  of  the  pure- 
state  Wrepresentability-problem,  where  one  has  to  de¬ 
cide  whether  the  reduced  density  operators  arise  from 
a  pure  state  of  N  fermions?  In  that  case,  the  ver¬ 
ifier  must  be  able  to  convince  himself  that  the  state 
he  gets  is  pure.  This  can  be  done  when  he  gets  two 
states  p  and  a  that  are  promised  to  be  uncorrelated, 
i.e.  that  he  gets  the  state  p  0  er:  then  Arthur  can  cal¬ 
culate  the  expectation  value  of  the  observable  tr(pcr), 
and  this  can  only  be  close  to  1  when  p  is  pure  and 
a  —  p.  Indeed,  if  tr(cr2)  <  1  —  e,  then  for  all  states 
r,  tr(crr)  <  ^/tr(tr2)  tr(r2)  <  (1  —  e)1/2  <  1  —  e/2.  The 
problem  however  is  that  Merlin  can  cheat,  and  hand  out 
a  correlated  state  to  Arthur,  such  that  the  above  test 
becomes  inconclusive.  This  is  precisely  the  feature  that 
distinguishes  the  complexity  class  QMA(k)  from  QMA 
f3:  in  QMA(k),  the  verifier  is  promised  to  get  a  tensor 
product  of  different  states,  and  it  has  been  conjectured 
that  QMA(k)  is  strictly  larger  than  QMA.  The  above 
discussion  shows  that  the  pure-state  ./V-representability 
problem  is  contained  in  QMA(k),  and  it  is  hence  plausi¬ 
ble  that  it  is  harder  than  the  ./V-representability  problem. 
It  would  be  very  interesting  to  investigate  whether  the 
problem  is  QMA(k)-complete. 

It  is  remarkable  that  checking  consistency  of  2-body 
reduced  density  operators  of  fermionic  states  is  so  hard, 
while  checking  consistency  of  1-body  reduced  density  op¬ 
erators  is  simple  @ .  This  can  be  easily  understood  from 
the  previous  discussion:  the  extreme  points  of  the  con¬ 
vex  set  of  1-body  density  operators  ( a\aj )  correspond 
to  ground  states  of  Hamiltonians  only  containing  bilin¬ 
ear  terms  in  a|  and  ay,  such  Hamiltonians  can  easily  be 
diagonalized  as  they  represent  systems  of  free  fermions, 
and  hence  the  consistency  problem  can  easily  be  solved. 
As  shown  in  2|,  consistency  can  be  decided  in  that  case 
based  solely  on  the  eigenvalues  of  the  reduced  density 
operators.  A  number  of  related  problems,  where  con¬ 
sistency  only  depends  on  the  eigenvalues,  have  been  in¬ 
vestigated  recently  13;  most  remarkably,  a  characteri- 
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sation  has  been  obtained  for  the  polytope  of  1-particle 
marginals  that  are  IV-representable  under  the  condition 
that  the  iV-particle  wavefunction  must  be  pure  0- 
These  results  have  to  be  contrasted  with  our  problem 
of  finding  the  IV-representable  2-body  density  operators, 
where  the  eigenvalues  alone  are  not  enough  to  decide 
consistency  but  also  the  eigenvectors  are  relevant.  Ac¬ 
tually,  let  us  consider  the  simpler  problem  of  deciding 
IV-representability  of  2-fermion  density  operators  where 
only  the  diagonal  elements  =  (aJatajCij}  are  speci¬ 
fied.  If  we  consider  the  case  d  =  2N  and  the  mapping 
discussed  above,  one  easily  finds  that  the  extreme  points 
of  this  set  would  be  obtained  by  ground  states  of  local 
spin  Hamiltonians  which  only  contain  commuting  az  op¬ 
erators.  These  correspond  to  spin-glasses,  and  so  the 
problem  of  deciding  IV-representability  of  {dy}  is  NP- 
hard  0  It  was  indeed  pointed  out  a  long  time  ago  that 
IV-representability  restricted  to  the  diagonal  elements  is 
equivalent  to  a  combinatorial  problem  jl9j|  that  was  later 
shown  to  be  equivalent  to  the  NP-hard  problem  of  decid¬ 
ing  membership  in  the  boolean  quadric  polytope  0- 
Let’s  finally  discuss  the  relevance  of  the  above  results  in 
the  context  of  quantum  chemistry.  We  have  shown  that 
it  is  a  hopeless  task  to  determine  the  ground  states  of 
all  local  fermionic  Hamiltonians,  and  in  particular,  that 
an  approach  by  means  of  the  IV-representability  prob¬ 
lem  is  intractable,  even  on  a  quantum  computer.  But 
it  is  possible  that  other  physical  systems,  e.g.  systems 
with  different  particle  statistics,  additional  symmetry  or 
a  limited  number  of  modes,  might  allow  for  an  efficient 
characterisation  of  the  two-particle  reduced  density  ma¬ 
trices,  and  hence  an  efficient  calculation  of  the  ground 
states  of  local  Hamiltonians.  This  seems  to  be  the  case 
in  e.g.  one-dimensional  translational  invariant  spin  sys¬ 
tems,  where  the  density  matrix  renormalization  group 
(2l|  allows  for  a  systematic  approximation  of  the  allowed 
convex  set  of  reduced  density  operators  from  within  |22l| . 
A  very  different  numerical  approach  has  been  developed 
in  the  context  of  quantum  chemistry  and  it  known  as  the 
the  contracted  Schrodinger  equation  Q.  In  that  method, 
one  approximates  the  convex  set  from  the  outside,  and 
the  IV-representability  problem  is  a  crucial  ingredient  of 
the  algorithm.  The  foregoing  discussion  shows  that  this 
approach  has  to  break  down  in  the  most  general  case,  and 
it  would  be  very  interesting  to  investigate  the  conditions 
under  which  those  approximations  are  justified. 

In  conclusion,  we  investigated  the  problem  of  IV- 
representability,  and  characterized  its  computational 
complexity  by  showing  that  it  is  QMA-complete.  Ob¬ 
viously,  the  theory  of  quantum  computing  was  a  prereq¬ 
uisite  to  pinpoint  the  computational  complexity  of  this 
classic  problem.  . 
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[24]  This  includes  search  and  optimization  problems,  e.g., 
“find  the  ground  state  energy  of  H.”  These  problems 
may  be  QMA-hard,  even  though  they  are  not  phrased 
as  yes-or-no  questions. 

[25]  Another  option  would  be  to  represent  each  qubit  %  us¬ 
ing  either  zero  or  two  fermions,  e.g.,  for  a  single  qubit, 
|0)  and  |1)  correspond  to  |fl)  and  aj&j|fl).  Define  the 

mapping  of  ^aj&;  +  bja;)  +  (biCii  +  aj&jj  ,  <jV  <-> 

i  (b\cii  —  affeij  +  i  (a\b\  -  fc;a;j ,  of  1  —  2 b\bi.  Then 
the  Hamiltonian  would  act  identically  on  the  subspace 
with  one  fermion  per  site,  and  on  the  subspace  with  zero 
or  two  fermions  per  site,  and  even  a  small  constant  in 
front  of  the  projectors  P;  would  be  enough  to  guarantee 
that  we  end  up  in  the  right  subspace. 

[26]  One  possible  set  of  observables  is  the  following:  First, 
define  a/  =  a;2aq,  for  all  pairs  I  =  {11,12},  *i  <  *2- 
Also  fix  an  ordering  on  the  pairs  I.  We  now  define  the 
following  observables:  Xu  =  a\aj  +  aja/,  for  all  I  -<  J; 
Yu  =  —ia\aj  +  ia'jai,  for  all  I  -<  J;  and  Zi  =  aja/,  for 
all  I  except  the  last  one.  These  operators  are  Hermitian, 
with  eigenvalues  in  the  interval  [—1,  1].  Taking  real  linear 
combinations,  these  operators  form  a  basis  for  the  space 
of  2-fermion  density  matrices. 

[27]  It  is  easy  to  see  that  K  is  contained  in  a  ball  of  radius 
R  =  VI,  since  for  all  a  £  K,  we  have  —  1  <  as  <  1. 

The  second  claim,  that  K  contains  a  ball  of  radius 
r  =  1/ poly(T),  is  less  trivial.  The  proof  is  as  follows. 

We  will  consider  A-representability  for  different  values 
of  N ;  let  Km  denote  the  set  of  all  vectors  a  that  are 
A-representable.  Obviously,  A'2  contains  a  ball  of  radius 
1/  poly(f')  (this  is  the  trivial  case).  We  will  prove  that  A'jv 
contains  a  ball  of  radius  1/  poly(£),  for  all  2  <  N  <  d  —  2. 

We  define  “particle-hole”  observables,  by  replacing  a; 
with  at,  and  vice  versa:  X\j  =  a/aj  +  aja\\  YdJ  = 
— ia/aj  +  iaja\;  and  Z\  =  aia\.  Let  a  denote  a  vector 
containing  expectation  values  for  these  observables;  let 
K'n  be  the  set  of  all  a  that  are  TV-representable. 

First,  we  claim  that  K2  =  K'd_2.  Notice  that  there  is  a 
natural  correspondence  between  the  2-particle  Slater  ba¬ 
sis  states  and  the  ( d  —  2)-particle  Slater  basis  states:  the 
2-particle  state  with  modes  i  and  j  occupied  corresponds 
to  the  (d  —  2)-particle  state  with  modes  i  and  j  empty. 

So  take  any  point  a  £  K2,  which  represents  the  ex¬ 
pectation  values  of  the  2-particle  observables  for  some  2- 
particle  state  a.  Use  a  to  construct  a  (d— 2)-particle  state 
r,  by  replacing  each  2-particle  Slater  basis  state  with  the 
corresponding  (d— 2)-particle  Slater  basis  state.  Then  the 
expectation  values  of  the  2-particle  observables  for  a  are 
exactly  the  expectation  values  of  the  2-hole  observables 
for  r.  So  a  is  in  K'd_2.  This  shows  that  A'2  C  K'd_2.  A 
similar  argument  shows  that  K'd_2  C  A2.  So  K2  =  Kd_2. 

Next,  we  claim  that  there  is  an  invertible  linear  trans¬ 
formation  A  that  maps  Kd-2  to  Kd-2. 

First  we  map  Kd_2  to  A'</_2.  Observe  that  we  can 
write  each  2-particle  operator  as  a  linear  combination 
of  2-hole  operators.  (This  holds  provided  we  restrict 


the  operators  to  act  only  on  the  subspace  of  (d  —  2)- 
particle  states.)  For  instance,  when  I  (~l  J  =  0,  a\aj  = 
aja\.  When  \I  n  J|  =  1,  we  write  equations  such  as 
a\a\aia,j  =  ajajaja;  =  —aja  J(1  —  aiaj)  =  —a,ja\  + 
a,ja\aia\  =  —a,ja\  +  aiajalaj;  then  use  the  fact  that 
aja\  =  (Ek&j  afc4)°j°i  =  Ek&j  a*ajaia l  When 
I  =  J,  we  have  a\ai  =  J}z,n/=0  aiai,-  (F°r  each  of  these 
equations,  it  is  easy  to  check  that  the  left  and  right  sides 
act  identically  011  all  (d  —  2)-particle  Slater  basis  states.) 

Thus,  for  any  (d  —  2)-particle  state  a,  the  expectation 
values  of  the  2-particle  observables  are  linear  functions  of 
the  expectation  values  of  the  2-hole  observables.  Thus  we 
have  a  linear  transformation  A  that  maps  Kd_2  to  Kd-2  ■ 

Similarly,  we  can  map  Kd-2  to  Kd_2.  We  write  each 
2-hole  operator  as  a  linear  combination  of  2-particle  oper¬ 
ators  (again,  restricting  the  operators  to  act  only  on  the 
subspace  of  (d— 2)-particle  states).  For  instance,  when  in 
J  =  0,  aia^j  =  a^a/.  When  |/n  J\  =  1,  we  write  equations 
such  as  aiaia^a]  =  ma^aia]  =  —  a]aj(l  —  a\ai)  =  — a]ai  + 
aja;a[ai  =  —  a^a;  +  a^ajaiai;  then  use  the  fact  that 
aja;  =  (^3  aj,afc)aja;  =  -d^Yhk^i,ja]a‘kakai- 

When  /  =  J,  we  write  a/a  j  =  aqaj  a»2at  =  (1  — 
a!1a*i)(l  ~  a j2 a;2 )  =  1  -  a^aq  -  aj2a;2  +  aja/;  then 
use  the  fact  that  1  =  ( d l]2)  X^_La/,°£:  alld  aja;  = 

— 3  ajajaja;  =  ajajaja;. 

Thus,  for  any  (d  —  2)-particle  state  a,  the  expectation 
values  of  the  2-hole  observables  are  linear  functions  of  the 
expectation  values  of  the  2-particle  observables.  Thus  we 
have  a  linear  transformation  B  that  maps  Kd-2  to  Kd_2, 
and  B  =  A-1. 

We  claim  that  Kd_2  contains  a  ball  of  radius  1/ poly(£). 
We  know  that  K'd_2  contains  a  ball  of  radius  l/poly(f?) 
(since  A'2  =  K'd_2),  and  we  will  show  that  the  map  A 
does  not  shrink  this  too  much.  Write  the  singular  value 
decomposition  A  =  UDV,  where  U  and  V  are  unitary, 
and  D  is  diagonal,  with  diagonal  entries  Du  >  0.  Let  B  = 
A-1.  Looking  at  the  matrix  elements  of  B,  we  can  see 
that  tr(BtB)  =  JW  |-Bqj2  <  poly(£).  At  the  same  time, 
tr(BfB)  =  tn{UD-1VV-1D-1U-1)  =  tr(AT2)  >  Djj2, 
for  all  i.  So  we  have  Du  >1/ poly(^),  for  all  i.  Thus  the 
map  A  shrinks  the  radius  of  the  ball  by  at  most  a  poly(f?) 
factor.  So  Kd-2  contains  a  ball  of  radius  1/ poly(^). 

Now  we  are  almost  done.  Notice  that  Ajv-i  3  Kn, 
for  all  3  <  N  <  d  —  2  (since,  if  a  is  consistent  with 
some  IV-fermion  state  a,  then  it  is  consistent  with  the 
( N  —  l)-fermion  state  a'  that  results  from  tracing  out 
the  last  particle).  So  we  conclude  that  A'jv  contains  a 
ball  of  radius  1/ poly(^),  for  all  2  <  N  <  d  —  2.  This 
completes  the  proof. 

[28]  For  example,  let  A;  and  \0i )  be  the  eigenvalues  and  eigen¬ 
vectors  of  O,  and  note  that  they  are  easy  to  compute. 
Add  an  ancilla  qubit,  perform  the  unitary  operation 
U  :  \0j)  |0)  1 — *  1 6i)  (Vl  —  A; |0)  +  x/A7|  1)) ,  and  measure 
the  ancilla  in  the  0,1  basis. 


